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Finite-dimensional wave turbulence refers to the chaotic dynamics of interacting wave 'clusters' 
consisting of finite number of connected wave triads with exact three-wave resonances. We examine 
this phenomenon using the example of atmospheric planetary (Rossby) waves. It is shown that 
the dynamics of the clusters is determined by the types of connections between neighboring triads 
within a cluster; these correspond to substantially different scenarios of energy flux between different 
triads. All the possible cases of the energy cascade termination are classified. Free and forced chaotic 
dynamics in the clusters are investigated: due to the huge fluctuations of the energy exchange 
between resonant triads these two types of evolution have a lot in common. It is confirmed that 
finite- dimensional wave turbulence in finite wave systems is fundamentally different from kinetic 
wave turbulence in infinite systems; the latter is described by wave kinetic equations that account 
for interactions with overlapping quasi-resonances of finite amplitude waves. The present results 
are directly applicable to finite-dimensional wave turbulence in any wave system in finite domains 
with 3-mode interactions as encountered in hydrodynamics, astronomy, plasma physics, chemistry, 
medicine, etc. 
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I. INTRODUCTION 

A. Weak-wave turbulence in finite-size systems 

" Wave turbulence" refers to the chaotic dynamics of 
nonlincarly coupled oscillatory modes. The phenomenon 
appears in a variety of physical context from surface wa- 
ter waves, through atmospheric planetary waves, plasma 
waves, acoustic waves in solids and fluids etc. Depend- 
ing on the strength of nonlinear interaction one distin- 
guishes weak-wave turbulence from strong-wave turbu- 
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lence. Weak-wave turbulence is characterized by a small- 
ness parameter £ which is roughly the root mean square 
of the ratio of the nonlinear to the linear term in the equa- 
tion of motion. For surface waves £ is about the ratio of 
the wave amplitude to the wave-length A, for sound in 
continues media this is the ratio of the density variations 
to the mean density, etc. 

The theory of weak- wave turbulence is particularly well 
developed in the limit of infinite systems where the ra- 
tio of the system size L to the characteristic wave length 
A is very large, L/X —> oo. In that limit the observed 
energy spectrum (energy distribution between modes) is 
well described by the so called "wave kinetic equations" 
that received considerable attention in the last half cen- 
tury, see e.g. [1—5]. We will call this regime of weak 
wave turbulence "kinetic wave turbulence" to distinguish 
it from other regimes, 11 finite- dimensional wave turbu- 
lence" and "mesoscopic wave turbulence" that will be 
introduced later. 

Notice that when the parameter L/X is of the order of 
unity the dynamics of waves can be very well described 
by low-dimensional chaotic models of the type studied in- 
tensively in recent decades , see e.g. [6, 7]. In this paper 
we explore the non-linear dynamics of weakly interacting 
waves when the parameter L/X is neither of order unity 
nor very large. This regime of parameters cannot be de- 
scribed either by kinetic equations or as low-dimensional 
chaos; it calls for new approaches and novel concepts, as 
partially demonstrated in this paper. 

To clarify the possible new regimes of weak-wave tur- 
bulence in various finite-domain systems we consider the 
general mathematical framework that takes the form of 
an energy conserving partial differential nonlinear equa- 
tion for a field ^(r, t). In a finite domain S one expands 
ty(r,t) in a complete set of eigenfunctions §>j(r) of the 
linearized dynamics that satisfy the boundary condition 
on the boundary dS: 

nr,t)=J2Mt)^(r) , (1) 

i 

where in general j can be a multiple index and the am- 
plitudes Aj(t) are functions of time but not of space. 
Accordingly the dynamics can be represented by a set 
of ordinary differential equations for the vector of ampli- 
tudes A(t) = {Aj(t)}, of the form 

^&=iu j A j (t)+NL j {A) (2) 

where ujj is the (real) eigenfrcquency of the jth mode; 
The symbolic term NL stands for the nonlinear contri- 
butions in this equation. According to the Poincare and 
Poincare-Dulac theorems [8] the nonlinear contributions 
can be brought to a normal form by a nonlinear change 
of variables. The nonlinear monomials that survive the 
change of variable are the resonant ones. The n-tuple 
(loi, . . .ojn) of cigenfrcqucncies is said to be resonant if 
there exists a relation of the form 

ujj = miOJi + to 2 o>2 + • • ■ rn n oj n . (3) 



The order of the resonance is X)fc m fc- The resonant 
monomials are of the form A™ 1 . . . A™ n . 

For the case of weak nonlinearities we invoke the small- 
ness parameter £ <C 1 to discard all the higher order reso- 
nances, keeping only the lowest available order. In other 
words, if there exist solutions to the equation 

LUj = LU m + LU n , (4a) 

we keep only the resulting quadratic monomials also 
known as "three-wave interactions", which satisfy 
the conservation law (4a). For example, in space- 
homogeneous, scale- invariant, isotropic, infinite wave- 
systems, in which the dependence of the wave frequency 
u>(k) on the wave vector k = |fe|, ui(k) cx k a , the three 
wave resonances 

w(fci) +w(fc 2 ) =w(|fci + fc 2 |), (4b) 

are allowed if a > 1 [1]. For a < 1 Eq. (4b) has no 
solutions and one needs to account for higher order res- 
onances. 

In this paper we focus on problems for which Eq. (4a) 
has solutions, determining the leading nonlinearity. We 
note that as the ratio L/X increases, there may be more 
and more eigenfrcquencies that satisfy equation (4a). In 
particular, while at small values of L/X we can expect 
only isolated resonant triads of waves, for larger values 
of L/X triads can share a common mode and the num- 
ber of coupled resonant triads increases considerably, fi- 
nally forming infinite clusters of connected triads. Anal- 
ysis of the ensuing dynamics under the influence of such 
growing clusters is the main subject of the present pa- 
per. We will focus here on the case of small enough 
nonlinearity parameter £ to ensure that only waves with 
exact resonances are important. In this case one can 
consider only clusters of connected resonant triads of in- 
teracting waves. We will refer to the chaotic dynamics of 
interacting waves in this regime as "finite-dimensional 
wave turbulence" , to stress the importance of the fi- 
nite number of interacting modes with exact wave res- 
onances. With increasing of wave amplitudes one has to 
account also for quasi-resonanccs. This type of wave tur- 
bulence was called "discrete wave turbulence" [9] . In con- 
trast, in infinite systems, the resonance conditions [e.g. 
Eq. (4b)] has infinitely many solutions; then usually the 
(kinetic) wave turbulence can be described by wave ki- 
netic equations. A more detailed analysis [10] shows that 
in the plane (L/X , () there exists a region of parameters 
where there exists weak-wave turbulence whose proper- 
ties are intermediate between finite-dimensional and ki- 
netic regimes. Some features of this type of turbulence, 
called "mesoscopic wave-turbulence" , were observed, for 
example, in [4, 11]. 

To study finite-dimensional wave turbulence we focus 
here for concreteness on the example of the barotropic 
vorticity equation on a sphere; this is an idealized model 
for atmospheric planetary (Rossby) waves [12], shortly 
described in Sec. II. Planetary-scale motions in the ocean 
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and atmosphere are due to the shape and rotation of the 
Earth, and play a crucial role in weather and climate 
predictability [12] . Oceanic planetary waves influence the 
general large-scale ocean circulation, can intensify the 
currents such as the Gulf Stream, as well as push them 
off their usual course. For example, a planetary wave 
can push the Kuroshio Current northwards and affect the 
weather in North America [13]. Atmospheric planetary 
waves detach the masses of cold or warm air that become 
cyclones and anticyclones and are responsible for day-to- 
day weather patterns at mid-latitudes [14]. 

Recently a new model [6] was developed for the intra- 
seasonal oscillations in the Earth atmosphere, in terms 
of triads of planetary waves whose eigen- frequencies solve 
Eq. (4a). The study of the complete cluster structure in 
various spectral domains shows that both for atmospheric 
[15] and oceanic [16] planetary waves indeed the size of 
the clusters increases with the growth of the spectral do- 
main. In large clusters one finds both large and small 
wavenumbers, meaning that the energy flux between very 
different scales becomes possible, bringing with it the 
hallmark of turbulence. Nevertheless, both numerical 
simulations [17, 18] and laboratory experiments [19] in- 
dicate that the dynamics of wave systems with interme- 
diate value of Lj A do not obey the statistical description 
provided by wave kinetic equations. This domain calls 
for a specialized investigation which is initiated in this 
paper. 



waves in a large spectral domain i, \m\ < 1000, presenting 
in Sec. IV A the total number of clusters consisting of 
one, two, three, etc. triads with different topologies. The 
histogram of all cluster distributions with respect to the 
triad number in clusters, ranging from 1 to 3691 is also 
presented. In Sec. IV B we use a notion of PP-irreduciblc 
clusters [7] (important for the discussion of the energy 
flux in finite-dimensional wave turbulence) and show the 
histogram of their distribution in size, ranging from 1 to 
130 triads in the PP-irreducible clusters. 



Sections V and VI are devoted to numerical simulation 
and preliminary analytical studies of finite-dimensional 
wave turbulence in clusters typical for numerous physi- 
cal systems including the planetary waves described in 
Sees. Ill and IV. In Sec. V we begin with the anal- 
ysis of free evolution in small clusters: butterflies in 
Sec. VB and triple-triad clusters (stars and triple chains) 
in Sec. VC. The main questions, that we discuss in this 
Section are: 

- reasonable choice of initial conditions, interaction coef- 
ficients and data representation that allows to shed light 
on the typical features of finite-dimensional wave turbu- 
lence, that depend on many parameters; 

how the energy flux between traids depends on the 
type of connections, on the interaction coefficients, on 
the initial conditions and on the cluster topology. 



B. Structure of this paper 

Section II reviews the properties of atmospheric plan- 
etary waves which are important for our analyis: in 
Sec. II A we consider barotropic vorticity Eq. (5) on a 
sphere and its dynamical invariants (6); in Sec. II B wc 
project Eqs. (5) and (6) on the spherical basis; and in 
Sec. II C we analyze the properties of the resulting inter- 
action coefficients. 

In Sections III and IV we study the topology and 
other properties of finite size clusters of resonant tri- 
ads of planetary waves that influence the dynamics of 
finite-dimensional wave turbulence. In Sec. Ill we begin 
with small clusters of resonant triads. In Sec. Ill A we 
overview equation of motion (14) for isolated resonant 
triad and its dynamical invariants (15), present them in 
the Hamiltonian form (20), and use the notion of "ac- 
tive" and "passive" modes [7] denoted as A- and P- 
modes; these notions are crucial for our theory. In the 
next Sec. Ill B we move on to double-triad clusters, PP- 
, PA- and AA-butterflies, their Hamiltonian equation of 
motion (25), Hamiltonian (26) and Manlcy-Rowe dynam- 
ical invariants (27). Similar discussion for the triple-triad 
clusters, stars, chains and triangles, is given in Sec. Ill C. 
Specific for atmospheric planetary waves, 6-triad cluster, 
caterpillar, is presented in Sec. HID as an example of a 
more complicated cluster structure. 

In Sec. IV we study clusters of atmospheric planetary 



In Sec. VI wc study finite-dimensional wave turbulence 
with a constant energy flux in the long-chain clusters, 
consisting of a large number of triads (N S> 1). For this 
goal we introduce pumping of energy into the the leading 
(first) triad and damping in the driven (last) triad. We 
discuss in Sec. VIA how to mimic the energy pumping 
and energy damping in our particular problem and what 
are the necessary conditions of stationarity, Sec. VI B. In 
Sec. VI C we show that the distribution of mode ampli- 
tudes in long chains is universal in the following sense: 
it is asymptotically independent of number of triads in 
the chain in the limit of large N and of interaction coef- 
ficients (in a wide region of their definition). Moreover, 
the distribution for the forced case practically coincides 
with that for free evolution from initial conditions, cor- 
responding to the forced stationary case. 



In Sec. VII A we briefly summarize the main features 
of finite-dimensional wave turbulence discovered in this 
paper and formulate in Sec. VII B some important ques- 
tions in this field that remain unstudied. Our feeling 
is that the present paper presents many more questions 
than answers, and all these questions (and many other 
related ones) belong to a new field of study of weak-wave 
turbulence: finite-dimensional and mcsoscopic wave tur- 
bulence in finite-size physical systems. 
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FIG. 1: Color online. Resonant triad Ai, see Tab. I below. 
Wave numbers of modes m , £ are shown inside ovals. Red 
arrows are coming from an active mode A (with frequency 
ujs = uji + L02) and show directions of the energy flux to the 
passive P-modes (with frequencies uji and 0*2). 



II. ATMOSPHERIC PLANETARY WAVES 
A. Barotropic vorticity equation on a sphere 

Planetary waves in the atmosphere pose a rich and 
complicated problem which is influenced by the earth 
topography, the vertical temperature profiles (varying 
between land and ocean), global winds etc. We do 
not attempt here to take into account all this richness. 
The essence of the interesting dynamics can be gleaned 
from simplified models. A very simplified model of at- 
mospheric planetary waves (discussed, e.g. by Silber- 
man [20] ) is provided by the barotropic vorticity equation 
on a rotating sphere for the dimensionless stream func- 
tion ip(6, if, t). The variables t, and <p are the time, the 
latitude (-vr/2 <9<n/2) and longitude (0 < <p < 2tt) 
on the sphere. The equation reads [20]: 



dip 
A 



1 {dip d 
sk[0\dpd0 



dip 



1 



0_ 

86 



? V2ficos0 + AV> 
1 2 



sin 



1—) 

86 J sin 6* dtp 2 



(5) 



where A is the angular part of the spherical Laplacian 
operator. The stream-function gives rise to the velocity 
v = fl R [z x V^>] , where and R are the angular velocity 
and radius of the Earth and z is the vertical unit vector. 

Equation (5) conserves the energy E and the enstrophy 
H, which are defined by: 

E = i J dip J I W 1 2 sin 0d0, (6a) 


2-7T 7T 

H = dip J I A V> 1 2 sin 9d6 . (6b) 


B. Projection on the spherical basis 

The eigenfunctions of the linear part of Eq. (5) are 

*j = ^ (0, tp, t) = Yj (0, V ) cxp( i Uj t) , (7a) 



where the frequencies of planetary waves are: 



u{l h mj) = -2 mjQ, j £j(£j + 1) 



(7b) 



From now on we use the shorthand notation j = (£j, rrij) 
for the eigen- numbers £j , rrij of the spherical harmonic 
Y^{0,p) 

Yj = Y™ 3 {0, ip) = P™ 3 (cos 0) exp(i mj <p) , (7c) 

with the associated Lcgendrc polynomials Pj = 
PT 3 (cos 0), normalized as follows: 



Ppppsin0d0 = 5{t,t') 



P7 m = PI 



(7d) 



Here 6(£, £') is the Kronecker symbol (1 for £ = £' and zero 
otherwise). The integer indices m and (£—m) are the lon- 
gitudinal and latitudinal wave-numbers of the £, m-mode; 
they count the number of zeros of the spherical function 
along the longitudinal and the latitudinal directions. Be- 
low we refer to the range of m and £ as the "spectral do- 
main" . For the approximation of two-dimensional atmo- 
sphere to hold, the wave-length is supposed to be much 
smaller than the atmosphere's hight. If wc estimate the 
wave-length as the distance between the appropriate ze- 
roes of spherical function, the length of the equator at 
about 40 000 km and the height of the atmosphere at 
about 40 km, we understand that the approximation of 
a two-dimensional atmosphere holds up to £ < 1000. 

Expanding the function ip(0,ip,t) in the basis (7a) we 
get 



iP(6, i p,t)=J2 A j*j 



(7e) 



Substituting in Eq. (5) one obtains the governing equa- 
tions for the "slow" amplitudes Aj = A™ 3 (t) of the plan- 
etary waves: 



dAj l ■sr^ 

~dT = ^W : 2^ N ^ sZ 'j\r,sArA s 



2N. 



(8a) 



x exp[i(uj r + uj s — u)j)t]S(mj, m r + m s ) , 
Nj = £j{£j + l), N rtS = N r -N 8 , (8b) 



where the three-wave interaction coefficients are 



ryfrij \m r ,m s 



7T 

J Zj lr , s (e)de, 



(9a) 



with the "interaction integrand" 

dP? 



m r Pp 



m s P e 



dPl n 



do 



(9b) 



We note that Eq. (8) is an exact consequence of the 
barotropic vorticity equation, without any assumption 
about the existence of a small parameter. Among the 
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interactions appearing in this equations there are many 
non-resonant ones, in which the exponent exp[i(aj r + 
lu s — u>j)t] is not unity. All these interactions can be 
removed by a change of variables, which however will 
result in new nonlinear terms (higher than quadratic), 
see e.g. Sec. 1.1.4 in [1]. Assuming that the wave am- 
plitudes are small enough, all these can be disregarded, 
bringing the final equations back to the same form as 
in Eq. (8), but including only resonant triads for which 
cxp[i(aj r + uj s — Wj)t] = 1. 



III. SMALL CLUSTERS OF RESONANT 
TRIADS 

In this Section we formulate equations of motion and 
motion invariants of small clusters of resonant triads, 
study their topology and other properties, that affect 
on the dynamics of finite-dimensional wave turbulence. 
Brief analysis of these questions was given in Ref. [7]. 

A. Active and passive modes in a resonant triad 



C. Necessary conditions for non- vanishing 
interaction 



The first necessary condition that guarantees finite in- 
teraction amplitudes follows from the axial symmetry of 
the problem and is reflected in the Kronecker symbol in 
Eq. (8a): 



(10a) 



Second, the spherical symmetry of the nonlinear term in 
Eq. (5) leads to the conservation of the square of the 
total angular momentum of the system. This translates 
to the triangle inequality for vectors £j = £ r + £ s in each 
triangle with non-zero interaction amplitude: 

<£j <£ r +L ■ (10b) 
Third, the explicit form of Eq. (9a) requires that 



£j + i r + i s is odd. 



(10c) 



Otherwise integrand (9b) is odd function of cos 8 and 
integral (9a) is zero. 

It can be shown by integrating Eq. (9) by parts that 
whenever Eq. (10a) is fulfilled the three interaction co- 
efficients satisfy: 



mj \m r ,m a gm a \— m r ,mj 7 m r \irij, — m s 



(11) 



From this follows that the interaction coefficients Z... sat- 
isfy two Jacoby identities 



J j\r,s 



Zr\s,j ^s\j,r 



(12) 



NjZj\r,s + N r Z r \ S) j + N s Z s \j r — . 
As a result, Eqs. (8) have two integrals of motion 



(13a) 
(13b) 



which are nothing else but the energy (6a) and enstro- 
phy (6b), presented in the basis (7). 



In this paper we refer to a "resonant triad" whenever 
we have three modes (j, r, s) whose frequencies satisfy the 
triad resonance condition ujj = lu t . + uj s . Accordingly, the 
equations for the slow amplitudes of modes in resonant 
triads are written, (after relabeling according to r — > 1, 
s — > 2 and j — > 3) as follows: 



N 2 ^ = iN ltS ZA$A 3 , 

A A* 

N3 ti = iN ^ ZA l A *2 ■ 



(14) 



E = -(JVi|Ai| 2 +JV 2 |A a | a + JV s |^3| 2 ) 



H = -(Nt\Atf + Ni\A^ + N£\A 3 \ 



In this case the conservation laws (13) take the form: 

(15a) 
) • (15b) 

Taking for concreteness an example for which 

h >t s >t 2 , (16) 

we make in Eqs. (14) a linear change of variables Bi = 
ctiAi such that 



■WN li2 N lt 3/y/N 2 N 3 , 



Oil 

a 3 = 



(17) 



This results in equations with real coefficients that in- 
volve only one interaction amplitude Z: 



dB 1 

~dT 

dB 2 

dt 
dB 3 

dt 



ZB 2 B 3 , 
ZB\B 3 , 
—ZB\B 2 



(18) 



This is a dynamical system corresponding to the simplest 
possible resonant cluster. Equations (18) are symmet- 
ric with respect to replacing two low-frequency modes 
\ <f> 2. The mode with highest frequency (which in this 
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FIG. 2: Color online. Examples of isolated butterflies. All the notations are as in Fig. 1. In particular: red arrows are coming 
from the active modes A and show directions of the energy flux to the passive P-modes. The letters "a" and "b" in square 
boxes denote triads and will be used as subscripts in the corresponding evolution equations for amplitudes and for integrals of 
motion. In studies of free evolution in Sec. V the initial energy is concentrated in the "leading" a-triad in "individual" modes 
with amplitudes So and Bo and then goes to the "driven" &-triad. 



paper will be always denoted by subscript "3 ") is spe- 
cial. When Eq. (16) does not hold one can find a similar 
change of variable for any relations between the magni- 
tudes of the three indices £j. 

The system (18) has two independent conservation 
laws (known as Manley-Rowe integrals) 

'ha = \B 2 \ 2 + \B 3 \ 2 = (EN 1 - H)N 23 /N 1 N 2 N 3 , 

/l3 = |i?l| 2 - ln - 12 



B 3 \ 2 = (EN 2 - H^/N^Ns , (19) 



hi 



I 



13 



\Bi\ 



I -B 2 



which are linear combinations of the energy E and en- 
strophy H defined by Eqs. (15). 

Obviously, Eq. (18) can be written in the Hamiltonian 
form: 



dt 



dB* 



with the interaction Hamiltonian 



-ffint — lZ{B\B 2 B 3 — BIB2B3) 



(20a) 



(20b) 



which is an additional integral of motion. In terms of old 
variables Aj 



rr . ^1,2^1,3^3,2 , „ . .* 

H\„t = iZ (AiAoAn 



A{A* 2 A 3 ). (20c) 



By direct calculation it is easy to check that Eq. (20c) is 
an integral of motion of the dynamical system Eq. (14). 

On the face of it Eqs. (18) involves six dynamical vari- 
ables: i.e. B's and their complex conjugates. In fact, 
using the standard representation of the complex ampli- 
tudes Bj in terms of real amplitudes Cj and phases df. 



Bj = Cj exp(i9j) , 



(21a) 



one recognizes that the right-hand-side (RHS) of Eq. (18) 
depends only on a single combination of phases in the 



triad which affects the dynamics. We refer to this com- 
bination as the triad phase: 



<P = Vl 



The triad phase appears in Eq. (20b) as follows: 
H int = -2Z\B 1 B 2 B 3 \ sin <p . 



(21b) 



(22) 



Thus we have a four dimensional phase space with three 
integrals of motion, resulting in a simple periodic trajec- 
tory for almost all conditions. For more details see [6]. 
Nevertheless even this simple dynamics offers the first 
opportunity to discuss the energy flow within a cluster of 
interacting modes. 

To this aim we discuss the evolution of the triad of 
amplitudes with special initial conditions, when only one 
mode is appreciably excited at zero time. If B\(t = 0) 3> 
B 2 (t = 0) and B x {t = 0) > B 3 (t = 0), then I 23 {t = 0) < 
Ii 3 (t = 0). The integrals of motion are independent of 
time, therefore /13 3> I 23 at all later times, and hence 
\Bi(t)\ 2 > |£ 2 (i)| 2 . Moreover, |Bi(i)| 2 > |^ 3 (£)| 2 at all 
times. Indeed, the assumption |i?i(i)| 2 < |i?3(£)| 2 yields 
ii3 ~ I23, which is not tenable. This means that the 
Wi-mode, being the only essentially exited one at t = 
cannot redistribute its energy to the other two modes 
in the triad. The same is true for the W2-mode. For 
this reason we refer to the lower frequency modes with 
frequencies u>\ < lo 3 and uj 2 < lo 3 "passive modes", or 
P-modes. 

On the other hand, the conservation laws (19) can- 
not restrict the growing of P-modes from initial condi- 
tions when only W3-mode is appreciably excited. In this 
case the P-mode amplitudes will grow exponentially [22]: 
|Bi(i)| , \B 2 (t)\ oc exp[\ZB 3 (t = 0)\t] until all the modes 
will have comparable magnitudes of their amplitudes. 
Therefore we refer to the uj 3 -mode as an "active mode" , 
or A-mode. An A-modc, being initially excited, is capa- 
ble of shearing its energy with two P-modes within the 
triad. 
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B. Double-triad clusters: Butterflies 

An arbitrary cluster in our wave system is a set of 
connected triads. Examples of the simplest clusters, con- 
sisting of two triads connected via one common mode are 
shown in Fig. 2. They will be referred to as butterflies. 
Note that in principle one can have two triads connected 
by two modes, but such a cluster does not satisfy all the 
conditions stated above. 

The dynamics of a cluster depends on the type of the 
mode which is common for the neighboring triads. Corre- 
spondingly we can distinguish three types of butterflies: 
PP-, AP- and AA-butterflies, In this Section we consider 
the equations of motion, the invariants and the restric- 
tions on dynamical behavior, that follow from the exis- 
tence of invariants for relatively small clusters consisting 
of two triads. In the following Sees, we consider these 
questions for clusters consisting of three and six triads. 

Butterflies, as shown in Fig. 2, consist of two triads a 
and 6, with wave amplitudes Bji a , Bj\ b , j = 1,2, 3, 
connected via one common mode. For PP-butterfly the 
common mode is passive in both triads, say 

B 1 \ a = B 1 \ b PP-butterfly; (23a) 

for an PA-butterfly the common mode is passive in a- 
triad and active in the second, 6-triad: 

B 1{a = B 3{b PA-butterfly; (23b) 

while for an AA-buttcrfly the common mode is active in 
both triads: 

B 3 \ a = B 3 \ b AA-butterfiy . (23c) 

The equations of motion for these systems follow from 
the Eqs. (8) under the condition of small nonlinearity and 
from the resonance conditions in both triads: 

Ul\ a +U2\a =U 3 \a, LU 1]b + UJ 2 \ b = LU 3 \ b , (24) 

with the obvious requirement that the frequencies of 
the common modes are the same. After a change of 
variables similar to Eqs. (17) and elimination of one 
common mode, the resulting equations for PP-butterfly 
{Bi\a = B 1[b ) are: 

{B\\ a = z aB^ a B z \ a + Z b B*^ b B z \ bl , 
B 2 \ a = z aB^ a B z \ a , B 2 \ b = Z b B^ a B 3 \ b , (25a) 
B 3 \ a = -Z a B 1 \ a B 2 \ a , B 3 \ b = —Z b B 1 \ a B 2 \ b ■ 

For PA-butterfly (with B^ a = B 3 \ b ) they are 

Bi\ b = Z b B*^ b B 3 \ b , B 3 \ a = —Z a B 3 \ b B 2 \ a , 
< B 2 \ b = Z b B^ b B 3 \ b , B 2 \ a = Z a B*\ b B 3 \ a , (25b) 
,B 3 \ b = -Z b B x \ b B 2 \ b + Z a B*\ a B 3 \ a , 



and for AA-butterfly (with B 3 \ a = B 3 \ b ): 

B\\ a = Z a B*\ a B 3 \ a , Bi\ b = +Z b B*^ b B 3 \ a , 
B 2 \ a = Z a B^ a B 3 \ a , B 2 \ b = Z b B*^ b B 3 \ a , (25c) 

B 3 \a = -ZaB^aB^a - Z b B x \ b B 2 \ b . 

All these equations can be obtained from the canonical 
equations of motion (20a) using the Hamiltonian 

H int = 21m ^Z a Bl\ a B 2 \ a B 3 \ a + Z b Bl\ b B 2 \ b B 3 \ b j , (26) 

in which the conditions (23) have to be fulfilled for each 
particular butterfly. 

In addition to the Hamiltonian (26) we have three more 
invariants of the Manley-Rowe type for each butterfly. 
For the PP-, PA- and AA-butterflies they respectively 
are 

^2,3|a = \B 2 a \ 2 + \B 3 \ a \ 2 , ^2,3|6 = \B 2 \ b \ 2 + \B 3 \ b \ 2 , 

I\a,b = \Bi a ? + \B 3]a { 2 + \B 3lb \ 2 , PP; (27a) 

h,2\b = 1^1 b\ 2 - \B 2 \ b \ 2 , I 2 ,3|a = l^2|a P + l^3|a| 2 7 

I\a.b = \Bi b \ 2 + |S 3 | h | 2 + |B 3 | a | 2 , PA; (27b) 

4,2|a = \B± a \ 2 - \B 2 \ a \ 2 , A,2|6 = \B\\ b \ 2 - \B 2 \ b \ 2 , 

I\a.b = \B la \ 2 + \B 1{b \ 2 + \B 3{a \ 2 , AA. (27c) 

The first two invariants for the PP butterfly, /2,31a an d 
-^2,3|6> do not involve the common mode Bi\ a = Bi\ b , 
and are similar to the invariant I 23 , Eq. (19), for an iso- 
lated triad. We can make the following observation: If 
at t = the amplitudes in one triad exceed substan- 
tially the two remaining amplitudes of the butterfly, that 
is \Bx\a\, \B 2 \ a \, \B 3 \ a \ > |B 2 |6|,|-B 3 | 6 |, this relation per- 
sists. In other words, in PP-butterfly when any of the 
two triads, a or 6, has initially very small amplitudes, 
it is unable to absorb the energy from the second triad 
during the nonlinear evolution. 

The invariants /1,21b an d /2,31a f°r the PA butterfly 
do not involve the common mode B 3 \ b = Bi\ a ; they 
are similar to the corresponding integrals I\ 2 and I 23l 
Eqs. (19), for an isolated triad. If at t = the 6-triad is 
excited much more than the a-triad (and thus I\ y2 \ b S> 
-^2,3|a) the smallness of the positively defined invariant 
J 2 ,3|a prevents the a-triad from absorbing energy from 6- 
triad during the time evolution. The situation is different 
when the a-triad is initially excited and /2.3|a 3> Ii.2\b- 
In this case the initial energy of a-triad can be easily 
shared with 6-triad. The smallness of I\, 2 \ b only requires 
that during evolution « l-E^lbl- Under this type of 

the initial conditions we will call the a-triad "leading" 
triad, while the 6-triad will be referred to as "driven" 
triad. 

Finally, the invariants for the AA butterfly, I\_ 2 \ a and 
^i.2|b7 do not involve the common mode i?3| a = Bi\ b and 
are similar to I\ 2 , Eqs. (19), for isolated triad. Simple 
analysis of these integrals of motion shows that energy, 
initially held in one of the triads will be shared between 
both triads dynamically. 
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One of tree AAA-stars 




One of five AAP-stars 





One of six APP-stars 



One of eleven PPP-stars 



FIG. 3: Color online. Examples of isolated triple stars in the spectral domain m < £ < 1000. The letters "a", "b" and "c" 
in square boxes denote triads and will be used as subscripts in the corresponding evolution equations for amplitudes and for 
integrals of motion. All other notations as in Fig. 1 and 2. 



The conclusion that we can draw from these examples 
is general: any triad which is connected to any given clus- 
ter of any size whatsoever where the connection occurs 
via its passive mode cannot absorb the energy from the 
cluster, if initially the triad is not excited. In contrast, 
a triad connected to a cluster of any given size via an 
active mode can freely adsorb energy from the cluster 
during the nonlinear evolution. 



C. Triple-triad clusters: stars, chains and triangles 

Triple triad clusters consist of three triads, denoted as 
a-, b- and c-triads with the mode amplitudes denoted as 



B, 



Bj\ b , Bj\ c 



j = 1, 2, 3. There are three topologically 
different types of triple-triad clusters: with one common 
mode - stars, shown in Figs. 3; with two common modes - 
chains, and with three common modes - triangles, these 
clusters are shown in Fig. 4. Having in mind different 
types of common modes one distinguishes 13 types of 
triple-triad clusters, including four stars (AAA-, AAP-, 
APP-, PPP-stars, Fig. 3), 7 types of 3-chain-, and two 
types of triangle-clusters, in Fig. 4. All motion equa- 
tions can be written in the canonical form (20a) with the 



Hamiltoman 



.Hint - 21m { z aB*^ a B*\ a B 3 \ a + Z b B^ b B*^ b B : 



3\b 



z cBl\ c Bl\ c B 5 \A 



(28) 



in which one has to equate amplitudes of common modes. 

a. Stars have one common mode in three triads. 
For example, taking B 3 \ a = B 3 \ b = B 3 \ c in Eq. (28) one 
gets from Eq. (20) equations of motion for AAA-stars: 



= Z a B 2 \ a B 3 \ a , 


B 2 \a 


= Z a Bi\ a B 3 \ a 


= Z b B 2 \bB 3 \ a , 


B 2]b 


= Z a Bx\ B 3 \ a 


= Z c B 2 \ c B 3 \ [a , 


B 2 \ c 


= Z c B x \ c B 3 \ a 



(29a) 



B 



3|a 



Z a B 1 \ a B 2 \a - Z h B x \ h B 2 \ h - Z C B 1 \ C B 2 \ C 



Taking B 3 i a = B^ = Bi\ c one gets for AAP-stars: 



B\\ a — Z a B 2 \ a B 3 \ a , B 2 \ a — Z a Bi\ a B 3 \ a , 
Bub = Z h B 2 \ h B 3 \ a , B 2 \ h = Z a B x \ h B 3 \ a , 



L|6 
3|c : 
3\a 



Z c B 3 \ a B 2 \ c , B 2 \ c — Z c B 3 \ a B 3 \ c , 
-Z a Bi\ a B 2 \ a — Zi>Bi\ b B 2 \ b + Z C B 2 \ C B 



(29b) 



3|c 
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AP-PA triple chain cluster 




One of five AP-PP chains 




One of six PA-PA chains 






One of six AA-PA chains 



One of fourteen AA-PP chains 



One of sixteen PA-PP chains 






one of eighteen PP-PP chains 



One of two AA-PA-PP triangles 



One of two PA-PP-PP triangles 



FIG. 4: Color online. In the spectral domain m < t < 1000 there are 66 isolated triple-chain clusters of seven types (examples 
are shown in panels A — J-) and four triangle clusters with (AA-PA-PP)- and (PA-PP-PP)-connections (examples are shown in 
panels Tt , X)). Common PP-modes are split, denoting difficulty in the energy exchange between corresponding triads. Dashed 
lines separate PP-irreducible sub-clusters, discussed in Sec. IV B. 



Similarly one gets motion equations for: 



does not involve the common mode. For example, for 



PPA-star: 

-S3 1 a = " 
#3|& = " 



B\\ a — B 1 

— Z a B\\ a B 2 \ a , 

— ZbBi\ a B 2 \b , 
ZcB^ c B\\ a , 



B 



3|c 



B2\a — Z a B x \ a B 3 \ a 

B 2 \b = ZaB\\ a Bj,\b , 



B 



2|< 



ZcBuaB^, 



(29c) 



.Bi\ a — Z a B*\ a B 3 \ a + Z b B*\ b B 3 \ b — Z c Bi\ c B 2 \ 
PPP-star: B 1]a = B 1]b = B 1{c , , 

B 3 \ a = ^Z a B 1 \ a B 2 \ a , B 2 \ a = Z a Bi\ a B 3 \ a , 
B 3 \b = -Z b B x \ a B 2 \ h , B 2 \ h = Z a B x \ a B 3 \ h , 
B 3 \ c = -Z c B x \ a B 2 \ c , B 2 \ c = Z c B x \ a B 3 \ c , 
B\\ a = Z a B*^ a B 3 \ a + Z b B*\ b B 3 \ b + Z C B 2 \ C B 3 \ C 



(29d) 



In addition to Hamiltonian, all triple-star clusters have 
four invariants of the Manlcy-Rowe type, three of them 



PPA-star: 



h,2\a = l-^llal 2 - \B 2 \ a \ 



Z 2,3|& = 
I\a,b,c 



\B 2 \ b \ 2 - 
\B lla \ 2 



IB 



3|fcl 



IB 



3|o 



J 2,3|c = 

|B 3 | b | 2 



|b 2]c | 2 

- |B 3 | C | : 



IB, 



(30) 



Integrals I 2i3 \b and I 2t3 \ c prevent b- and c-triads (con- 
nected via the P-mode) from adopting energy of initially 
excited a-triad. In cases when b- and/or c-triad are ini- 
tially exited, a-triad can freely share their energy via con- 
necting A-mode. 

b. Triple-chains have two common modes in two tri- 
ads. As we mentioned, there are 7 types of triple-chains, 
that differ in type of connections, see Fig. 4. Similarly 
to the triple-star clusters, one gets equation of motion 
for triple-chains from the canonical Eq. (20a) with the 
Hamiltonian (28), in which one has to equate two pairs of 
amplitudes of common modes. For example, for PA-PA 
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chain, 



D. Caterpillar: 6-triad cluster 



PA-PA-chain: Bu a = B, 



3\b ■ 



But = B 



3 c : 



B\\ a 


— Z a B^ a B 3 \ a - 


- Z b B x \ b B 2 \ b , 


B 2 \a 


= Z " B l\a B 3\a , 


B3\a = —Z a Bi\ a B 2 \a, , 


B 2 \ b 


= Z b B*^ b B 3 \ b , 


(31b) 


B\\b 


= Z b B*^ b B 3 \ b — 


Z a Bi\ c B 2 \ c , 


Aw 


= Z c B^ c Bi\ b , 


B 2 \c = Z c B^ c Bi\ b . 



(31a) The largest cluster found in the domain £ < 21 con- 

sists of six resonant triads An . . . Ai6 with three PP-, 
one AP- and one AA-connection, see Fig. 5. The equa- 
tion of motion for this cluster, (called "caterpillar") can 
be obtained from Hamiltonian, similar to Eq. (28), but 
consisting of six terms: 



Again, besides Hamiltonian, all triple-chain clusters have 
four invariants of the Manley-Rowe type, two of them 
do not involve the common mode. In particular, PA-PA 
chain, governed by Eq. (31a) has the following invariants: 



16 

-ffcat = 21m ^ Z n B t\ n B l\ n B 3\n i 
n=ll 



(33a) 



in which we have to equate amplitudes of common modes: 



— i?l|12 , i?3|12 — -Bl|13 , i?2|12 — 



B 



3 14 



B 



3|15 ■ 



-Sl|14 — #1 



(33b) 



PA- PA chain: 



I\b,c — \ B l\b\ 2 - \B 2 \ b \ 2 + \B- ■ 



B l\a\ 2 
2 \A 



z 2,3|a 



|S 2 | a | 2 



|B 3 | a | 2 , h,2\c 



\ B 2\a? 



l^llc' 2 



3 2|6l 



(31c) 



c. Triple-triangles have three common modes in 
three triads, see Fig. 4H and X. Correspondingly, equa- 
tions of motion for AA-PA-PP and PA-PP-PP triangles 
one gets from the canonical Eq. (20a) with the Hamil- 
tonian (28), in which one has to equate three pairs of 
amplitudes of common modes. In particular, for AA-PA- 
PP triangle one has: 



B\\ a = -Bi|& , B 3 l a 



B 3\c 7 B 



316 



Bile, (32a) 



B\\ a — 
B 2 \ a = 
B 3\a = 
B 2 \ b = 



Z a B 2 \a B 3\a + Z b B*^B 3 \ b , 
Z a B ^ a B 3 \ a , 

— Z aBi\ a B 2 \ a — Z c B 2 \ c B 3 \ b . 

z bB{\ a B 3 \ b , 



(32b) 



B 3 \h — —ZhB-i \„B 



B 



3 3|6 
2|c 



6- D l|a- D 2|6 
Z cB*\ b B 3\a ■ 



Zc B 2 \ c B3\ a 



Triangle clusters have three Manley-Rowe invariants. In 
particular, AA-PA-PP triangle, governed by Eq. (32b) 
has the following invariants: 



Due to these five connections one has only 3x6 — 5 = 13 
complex equations for remaining amplitudes: 

Bi\n = Zn-B|| 11 5 3 | 11 + Zi2-B2|i 2 -B3|i2 , 
B 2 |n = ZnB^^B^n , i? 3 |ii = — Z n Bi\ 11 B 2 \ii , 
^2|i2 = Z\ 2 B*^ Y1 B 3 \ X2 + Zi4_B*| 14 J5 3 | 14 , 
B3112 = — Z\ 2 Bi\ 12 B 2 \ 12 + Zi 3 B^ 13 B 3 \ 13 

Bi\i4 = Zi 4 i?2| 1 2-0 > 3|i4 + Zi 6 B^ 1Q B 3 \ 16 j 
-83114 = Z U B 2 \\2 B ^\U + Z W B 2 \\a B 3\l& ' 
-Bl|l5 = Zi5B*^ 15 B 3 \ u , -B 2 |15 = Z 15 B *\\§ B 3\1A , 

B 2 \ 16 = Zi 6 B^ 14 B 3 \ 16 , i? 3 |i6 = —Zi(,Bi\ u B 2 \iq 
Caterpillar has six Manley-Rowe invariants: 

2 



(33c) 



z 2,3|ll 



l^2|ll I 2 



^2,3|13 = I 13 | 2 

h, 2 \15 = l^l|X5 I 2 

-^2,3|i6 = l^iiel 2 

-'111,12,13 



|B 3 |11| 
|B 3 |i 3 | 2 , 
IB2I15I 2 , 



IB 



(33d) 



Bmil 2 



l. -'|14, 15, 16 



|Bi 



1 14 1 



3|16l 
|B 3 ]n 

IB1115I 



| 2 + |B 3 |i 2 | 2 

|B 3 |15| 2 



|B 3 | 13 | 2 , 
|B 3 ii 6 | 2 ■ 



IV. HOW CLUSTERS ARE ORGANIZED 



AA-PA-PP triangle 

/1 



I\a,c — |B 2 | a | 2 



B 3 |a| 2 



6,c 



B 2 |&| 2 - |B 3 |b| 2 



/|a.6 = |B 2 | 6 | 2 -|S 1 | Q | 2 , 

|B 2 | C | 2 , 



>2|c 



(32c) 



Notice, that triple-stars and triple-chains have 10 real 
variables (seven amplitudes and three triad phases) and 
five invariants (Hamiltonian and four Manlcy-Rowcs'), 
while triple-triangles have only 9 real variables (six am- 
plitudes and three triad phases) and four invariants 
(Hamiltonian and three Manley-Rowes'). Therefore, all 
triple-triad clusters have five-dimensional effective phase 
space. Recall, that butterflies have three-dimensional 
phase space. Moreover, one can prove that any n-triad 
cluster has (2n — l)-dimensional effective phase space. 



This Section is devoted to the analysis of the structure 
of clusters of resonant triads of atmospheric planetary 
waves (based on the data-set of the exact solutions of 
resonance conditions u>j = uj r +uj s with restrictions (10), 
provided by E. Kartashova [23]). This analysis is impor- 
tant for the study of finite-dimensional wave turbulence 
throughout the present paper. A preliminary study of 
this issue can be found in Ref. [7] . 



A. "Meteorologically significant" clusters and their 
extension in large spectral domain 

Dealing with Atmospheric waves one learns that the 
" meteorologically significant" wave numbers arc believed 
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FIG. 5: Color online. Triads belonging to caterpillar are drawn by bold lines. New, connected to them triads, appearing in 
spectral domain m,l < 1000 are drawn by thin lines. The PP-reduction of the extended caterpillar to three triads (An, Aie 
and An. 4, according to Tab. I), PA-butterfly 1x112,13, AA-butterfly 1x114,15, and AAA-star is shown by dashed lines. The rest 
of notations as in Figs. 1 and 4. 



Clust. 


Ni 


Modes [m,£] 


N 2 


Connecting triads 


Ai 


1 


[4,12] [5,14] [9,13] 






A 2 


2 


[3,14] [1,20] [4,15] 


2.1 


[4,15] [10,24] [14,20] 








2.2 


[1,20] [14,29] [15,28] 








2.3 


[1,20] [15,75] [16,56] 


A 3 


3 


[6,18] [7,20] [13,19] 


3.1 


[2,15] [5,24] [7,20] 


A 4 


4 


[1,14] [11,21] [12,20] 


4.1 


[1,14] [9,27] [10,24] 


1x15,6 


5 


[2,6] [3,8] [5,7] 


5.1 


[4,14] [9,27] [13,20] 




6 


[2,6] [4,14] [6,9] 








7 


[6,14] [2,20] [8,15] 


7.1 


[2,20] [11,44] [13,35] 


1x17,8 


8 


[3,6] [6,14] [9,9] 


7.2 


[2,20] [30,75] [32,56] 








7.3 


[32,56] [26,114] [58,69] 


1x19,10 


9 


[3,10] [5,21] [8,14] 








10 


[8,11] [5,21] [13,13] 








11 


[2,14] [17,20] [19,19] 


11.1 


[2,14] [18,27] [20,24] 




12 


[1,6] [2,14] [3,9] 


11.2 


[6,44] [14,21] [20,24] 




13 


[3,9] [8,20] [11,14] 


11.3 


[9,35] [11,20] [20,24] 


Klll-16 


14 


[1,6] [11,20] [12,15] 


11.4 


[3,20] [45,75] [48,56] 




15 


[9,14] [3,20] [12,15] 








16 


[2,7] [11,20] [13,14] 







TABLE I: The first 3 columns provide the following data in 
the domain m,£ < 21 : the cluster's form, the triad numbers 
and the modes within a cluster; in the last two - the numbers 
of additional connecting triads and their modes that enlarge 
the corresponding cluster when the spectral domain \m\, I < 
1000 is regarded. 



to be limited to £ < 21 [21]. Nevertheless, as explained 
above, the spectral domain for the approximate two- 
dimensional atmosphere extends to £ < 1000. Counting 
explicitly how many clusters we have in this spectral do- 
main we find that there exists altogether 1965 isolated 
triads and 424 clusters consisting from 2 to 3691 con- 



nected triads. Among them there are 234 butterflies, 95 
triple-triad clusters, etc. (cf. the histogram in Fig. 6.4). 
For clarity of presentation we did not display in this his- 
togram the largest 3691-cluster, which we refer to as the 
monster. 

It can be seen that about 82.2% of all clusters are pre- 
sented by isolated triads and their dynamics has been in- 
vestigated in [6] in all details; the main findings are that 
the energy oscillates between the three modes in the tri- 
ads, with a period of oscillation that is much larger then 
the wave period. This period is inversely proportional to 
the root-mean-square of the wave amplitude. 

234 clusters in the spectral domain (~ 10.5%) are the 
butterflies discussed above; these are further analyzed 
below. Among them there are 131 PP-, 69 AP- and 
35 AA-butterflies. The 95 triple-triad clusters include 
25 "triple-star" clusters with one triple connection, sec 
Fig. 3. This set includes 3 AAA-, 5 AAP-, 6 APP- and 
11 PPP-stars. There are also 66 chain clusters with two 
pair connections and 7 combinations of the connection 
types, shown in Fig. 4. We found also four (two pairs) 
triple-triad clusters with three pair connections, belong- 
ing to two different types, see Fig. 4. 

Similar classification can be performed for all the other 
clusters. For example, the monster includes one mode 
(218,545), participating in 10 triads, three modes, par- 
ticipating in 9 triads, 5 modes - in 8 triads, 23 - in 7, 
50 - in 6, 90 in 5, 236 - in 4, 550 - in 3 and 1428 modes 
- in 2 triads (butterflies). The analysis of their dynami- 
cal behavior depends crucially on the connection type as 
shown above and detailed below. 



B. PP-reduction of larger clusters 

Large clusters can be divided into " almost separated" 
subclusters connected by PP-connections (e.g. Figs. 5, 
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A 




The 3691-triad cluster (monster) is not shown 



B 




9 11 13 15 17 19 21 

Size of sub-cluster 



Sub-cl. size 


21 


21 


28 


30 


32 


33 


35 


37 


38 


11 


121 


127 


130 


Amount 


2 


3 


1 


1 


1 


2 


1 


1 


1 


1 


1 


1 


2 



Amount of sub-cluster of size > 21 shown in the table 



FIG. 6: Color online. Horizontal axes denote the number of 
triads in the cluster while vertical axes show the number of 
corresponding clusters (panel A) and PP-irreducible subclus- 
ters (panel £>) . 



7 and 8). If such a sub-cluster cannot be divided fur- 
ther into smaller clusters connected by PP-connection 
we refer to it PP-irrcducible cluster. For exam- 

ple, triangle cluster in Fig. 4 I can be PP-reduced to 
a triad and a PA-butterfly, while clusters in Fig. 8 A and 
B can be PP-reduced to two and three individual triads 
respectively. The cluster in Fig. 7 is PP-reduced into 
three triads and PA-butterfly, while extended caterpil- 
lar in Fig. 5 is PP-reduced into three triads, PA-, AA- 
butterfly and AAA-star. One can see that the 14-triad 
cluster, shown in Fig. 9 A can be PP-reduced into 8 triads 
and 6-triad cluster. The 16-triad cluster (Fig. 9 B) can be 
PP-reduced into 5 triads, AA- and PA-butterfly, AAP- 
star and 4-triad cluster, consisting of AAA-star with one 
AP-connected triad. 

In the spectral domain £, \m\ < 1000 the largest PP- 
irreduciblc sub-cluster belonging to the monster consists 
of 130 triads and is shown in Fig. 10. The statistics of 
PP-irrcducible sub-clusters are presented in Fig. 6 B. 



We learn from this analysis that many clusters cannot 
carry energy flux through PP connections; their dynam- 
ics is naturally reduced to the dynamics of PP-irreducible 
sub-clusters. Therefore it is sufficient to study carefully 
the dynamics of these PP-irreducible clusters to under- 
stand the properties of any cluster. 



V. NUMERICAL ANALYSIS OF FREE 
EVOLUTION OF TYPICAL SUB-CLUSTERS 

Our goal is to describe the energy flux through reso- 
nant triads in the regime of finite-dimensional wave tur- 
bulence. In the first Subsec. VB we consider free evolu- 
tion of the smallest clusters - butterflies - from asymmet- 
rical initial conditions, in which only one triad is excited 
to high amplitudes, exceeding by orders of magnitudes 
the initial amplitudes in the other triad. The questions 
are how the energy flux from the energetic "leading" triad 
depends on the type of connection, on the ratios of the 
interaction coefficients etc. Subsec. V C is devoted to 
the free evolution of the triple-triad clusters: stars and 
chains from initial conditions in which only the leading 
a-triad is substantially excited, the levels of excitation 
of the two other triads are much smaller. All these ex- 
amples, and PA-PA-.. PA chains, studied in the next Sec- 
tion, can serve as building blocks of bigger clusters and 
the knowledge about the energy flux through them allows 
to qualitatively predict efficiency of the energy transfer 
through bigger clusters. 



A. Methodology and numerical procedure 

The equations of motion for all clusters were prepared 
for numerical analysis using a specially designed algo- 
rithm that allowed an automatic implementation for any 
cluster, given the number of triads and the connectivity 
table. This served to avoid human errors in implementing 
large sets of equations. 

The equations of motion for the free evolution of but- 
terflies and triple clusters are stiff and were integrated 
using a multistep adaptive method based on numerical 
differentiation formulas [24]. 

For each system of equations the accuracy of integra- 
tion was controlled by testing the conservation of the 
relevant integrals of motion. The integration parameters 
were adjusted to keep the standard deviation cr(J) be- 
low a given threshold for the duration of the numerical 
runs. Since the main parameter that affects the accuracy 
in terms of the conservation of the integrals of motion is 
the maximal allowed time step, the preliminary calcula- 
tions were carried out with a requirement < 10~ 5 ; 
actually in most cases a(I) < 10 -8 was achieved. 

The evolution starting from several hundred initial 
conditions was analyzed and several representative con- 
ditions were chosen for the study of energy transfer in 
the clusters. All the conclusions regarding the discov- 
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FIG. 7: Color online. Triads A7 and As, that belong to PP-butterfly 1x17,8 are drawn by bold lines. New, connected to them 
triads A7.1, A7.2 and A7.3, appearing in spectral domain m,£ < 1000 (for numeration of triads see Tab. I), are drawn by 
thin lines. The PP-reduction of this cluster to three triads A7, As and A7.1 (according to the notation in Tab. I) and one 
AP-butterfiy is shown by dashed lines. The rest of notations as in Fig. 1. 




FIG. 8: Color online. Examples of PP-reduction of pair connection, panel A and of triple connection, panel B. Common (PP- 
and PPP-) modes are split stressing the difficulty of energy exchange between the corresponding triads. Dashed lines separate 
PP-irreducible sub-clusters, discussed in Sec. IV B. The rest of the notations are as in Fig. 1. 




FIG. 9: Large PP-reducible clusters. Common PP- and PPP-modes are shown by empty circles. The rest of modes are denoted 
by full (blue) circles. As in previous figures, outgoing (red) arrows indicate A-modes. 
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FIG. 10: Color online. Largest PP-irreducible cluster in the spectral domain I, \m\ < 1000. 



ered dependencies were verified by control calculations 
with stricter accuracy requirements. 

The equation of motion for forced chain clusters, stud- 
ied in the next Sec. VI, were integrated by both adaptive 
methods [24] and by 4th order constant time-step Rungc- 
Kutta for better control of accuracy. By construction of 
the model, these equations required accurate description 
of the last triad in the cluster to ensure proper energy 
dissipation. The integration parameters were adjusted 
to reproduce this fastest evolution, and therefore were 
automatically suitable for all other triads. We verified 
the convergence of the resulting statistics with respect to 
all relevant parameters. 



B. Free evolution in butterflies 

The simplest topology that allows consideration of the 
energy flux between resonant triads is the double-triad 
clusters - butterflies. 



1. Initial conditions, choice of the interaction coefficients 
and data representation 

In this Subsection we show that details of the time evo- 
lution in butterflies are very sensitive to the initial con- 
ditions, which define the values of the dynamical invari- 
ants. Therefore a reasonable choice of initial conditions, 
allowing to shed light on a "typical" time evolution in a 



relatively compact form is not obvious. At initial time 
t = we assign most of the energy to two individual (not 
common) modes of one (leading) triad, denoted below for 
concreteness as a-triad. The initial amplitudes of these 
two modes we denote as Bo and Bq (Fig. 2). To study the 
influence of the energy distribution between these modes 
we will use two types of initial conditions: 

Type I : B = 3.9 + 0.50 i, B = 3.7 + 0.93 i , (34a) 
Type II : B Q = 5.3 + 0.50 i, B o = 0.9 + 0.93 i . (34b) 

Both distributions (34) are complex and normalized such 
that |S | 2 + |i?o| 2 ~ 30. The difference between Eq. (34a) 
and (34b) is that in Eq. (34a) both amplitudes are simi- 
lar, while in Eq. (34b) they are quite different. 

For different butterflies we choose in the leading triad: 

PP-butterfly with S 1]a (0) = Bi| 6 (0) : 

B 2 \ a (0) = Bq, B 3{a (0) = B ; (35a) 
AA-butterfly with B 3 \ a = B 3 \ b : 

B 1{a (0) = B Q , B 2W (Q) = Bo ; (35b) 
PA-butterfly with B x \ a = B 3 \ b : 

B 2 |a(0) = B , B 3la (0) = B , (35c) 

as it is shown in Fig. 2. 

The initial conditions in the driven 6-triad we choose 
the same for all types of butterflies. They have much 
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smaller initial amplitudes, for example: 

' B llb (0) = Bi, , B h0 = C (0.05 + 0.02 i) , 
< B 2 |6(0) = £ 2i0 i ^2,o = C (0.02 + 0.05 ») , (36) 
k B 3 | 6 (0) = S 3 , , B 3 , = C (0.10- 0.02 i) . 

To study the dependence of the energy flow between tri- 
ads on the initial level of excitation of the driven triad 
we vary the energy content in the driven triad changing 
the coefficient C in Eq. (36) , taking in addition to C = 1 
also C = 0.1 and C = 0.01. In all the further simulations 
C = 1 if else is not mentioned. 

In this way the initial conditions for all butterflies are 
as similar as possible and we can study the difference in 
time evolutions, caused by different types of connections. 

Last but not least are the interaction coefficients. For 
concrcteness we chose interaction coefficients correspond- 
ing to A14 and Aie (Z±4 as 75, and Z\e ~ 15) as proto- 
types and use either {Z a = 75 , Z b = 15} or vise versa 
and sometimes {Z a = Z b = 15}. Since the change in 
the interaction coefficient re-normalizes the correspond- 
ing time scale, only their ratio is important for the dy- 
namics. In our case these ratios arc = 5 , 1/5 or 1; this 
allows us to study the butterfly dynamics with very dif- 
ferent values of the interaction amplitudes, which is typ- 
ically the case. Having in mind that special choices of 
the ratios of interaction coefficients may lead to integra- 
bility of clusters of resonant triads [25] (and see also [26] ) 
we verified that small variations of these ratios does not 
changed our conclusions concerning energy transfer in 
clusters. 

Recall that any butterfly has three quadratic integrals 
of motion, that involve only squares of five amplitudes of 
modes and therefore only 2 combinations of them are in- 
dependent. For the presentation we chose such combina- 
tions that are orthogonal to the corresponding invariants. 
Namely, for PP-butterffies, connected via B\\ a = B^b 
modes: 

•h,3\a = \B 2 \a\ 2 - l^3|a | 2 ! ^2.3|b = l^2|h| 2 ~ I ^3 1 fo 1 2 ! 

(37a) 

for AA-butterffies, connected via B 3 \ a = B^j, modes: 

■h.2\a = \B 1{a \ 2 + \B 2{a \ 2 , J 1)2 | b EE \B 1[a \ 2 + \B 2lb \ 2 , 

(37b) 

and for AP-butterflies, connected via B^ a = B 3 \ b modes, 

^2,3|a and Ji,2|fc- 

Time evolutions for these three types of butterflies with 
various initial conditions and choices of the ratio Z a /Z b 
(5 or 1/5) are shown in Figs. 11 - 15. Following Sub- 
section V B 3 is devoted to discussion of these numerical 
results. 

2. Effect of the type of connections and of the ratio Z a l Z b 

a. PP-butterfly has the most trivial time-evolution, 
see Fig. 11 for Z a = 75, Z b = 15, (panel A) and 



Z a = 15, Z b = 75, (panel B). As expected, there is 
practically no energy exchange between triads: ampli- 
tudes J 2 ,3\ a i ^2,3|&> defined by Eq. (37a), oscillate within 
the domains, that are determined by initial conditions. 
We show these evolutions just to illustrate our analytical 
result that PP-connection is non-penetrative for energy 
in both directions at any time. 

b. AA-butterfly is a promising candidate for the en- 
ergy transfer. Time evolutions of Ji i2 \ a (t) and Ji i2 \ a {t) 
for AA-butterfly [defined by Eqs. (37b)] are shown in 
Fig. 11 £ , T . Indeed, as one sees in panel T the peaks of 
the amplitudes J\ y2 \ a {t) in the leading triad and Ji :2 \b(t) 
in the driven triad are close to 30 for Z a = 15, Zb = 75. 
Quite unexpectedly, the energy transfer may be not effi- 
cient even with AA connections. Indeed, as one sees in 
Fig. 11 £, when the leading a-triad has larger Z a = 75, 
the sum |i?i| 2 + I-B2I 2 hi the leading triad oscillates be- 
tween as 30 and as 20 being close to its initial value as 30. 
At the same time, this sum in the driven triad oscillates 
close to its initial value as 0.01. Therefore the energy 
transfer is strongly suppressed if Z a Z b . 

To understand the difference between these two cases, 
consider time dependence of the common mode B 3 (t), for 
the case Z a = 15 , Z b = 75 shown in Fig. 12 A. One sees 
fast oscillations of the common mode B 3 (t) with (almost) 
zero mean and some frequency £l a that can be estimated 
as Z a B a , where B a = J\_ 2 \ a is a characteristic value 
of mode magnitudes in a— triad. An illustration of this 
estimate one sees in Fig. 11 £, where (for the same value 
of \B a \) the oscillation frequency of the triads with Z a as 
75 is much higher than that for triads with Z a as 15. 
As one can see from the equations of motion (25) for 
butterflies, the mean energy flux from the leading, a—, 
to the driven 6-triad (with A-conncction) , e a ^ bl can be 
written as 

e a ^b = 2Z b Rc[B 1 \ b B 2 \ b B^ b ] oc cos (j> b , (38) 

where ip b = (p±\ b + ip 2 \ b — ips\ b is the triad phase, intro- 
duces by Eq. (21b). Fast oscillations of B a (t), equivalent 
to fast growth of V 3 3|b(0 with the speed d<p 3 \ a (t) / dt = 
d<P3\b(t)/dt as n a (t) , lead to self-averaging of cos^b al- 
most to zero, if phases ipi\ b and tp 2 \ b cannot react fast 
enough to variations of ip 3 \ b {t). This is a qualitative ex- 
planation of the observed fact that the energy flux from 
a— to 6— triad is strongly suppressed if Z a 3> Z b and 
is observable (but still suppressed) if Z a < Z b . One can 
hope that the a— b energy flux can be significant if -B3| a (t) 
has a nonzero mean and therefore (<fi3\b(t)} also does not 
vanish. But this is impossible under the initial condi- 
tions of interest: we do not want to put initially energy 
to b— triad, taking large value of B 3 \ a (t) = B 3 \ b (t) at 
time zero. Therefore at t = | 1 a | |-Biu|, |-B 2 i a |. 
If so, according to the equation of motion, the time 
derivative B 3 \ a = —Z a B\\ a B 2 \ a is large enough to allow 
B 3 \ a to cross quickly zero and to reach significant value 
with different sign. Having in mind periodical character 
of evolution of isolated triad it practically means that 

(B 3 \ a ) « v^TF)- 



16 



A 



PP butterfly, leading triad, Z = 75 




1 2 3 4 5 

xKf 3 PP butterfly, driven triad, ^=15 




Blockade of the energy transfer in PP-butterfly 



-26 
J -28 
-30 

0.02 r 

-0.02 



B 



PP butterfly, leading triad, Z = 15 



1 2 3 4 5 

PP butterfly, driven triad, 7^= 75 

illlllilllllllii,. .iiiflllllllfiii,, „ 



3 
time 



Blockade of the energy transfer in PP-butterfly 



PA butterfly, leading triad, Z= 75 




10 



1 2 3 4 5 6 

* 10 PA butterfly, driven triad, Z= 15 




V 



PA butterfly, leading triad, Z=15 




1 2 3 4 5 

PA butterfly, driven triad, Z= 75 



3 
time 



llliiiil ' , Lli 

4 5 6 



Strong suppression of the energy transfer in PA-butterfly 



Efficient energy transfer in PA-butterfly 



40 r 



3 20 



0.015r 
£ 0.01 

r 

0.005 



£ 

AA butterfly, leading triad, Z= 75 




AA butterfly, leading triad, Z=15 



1 2 3 4 5 

AA butterfly, driven triad, Z= 15 



%W '"iw' %¥ l i 




1 2 3 4 5 

AA butterfly, driven triad, Z= 75 



3 
time 




Strong suppression of the energy transfer in AA-butterfly 



Efficient energy transfer in AA-butterfly 



FIG. 11: Color online. Time evolution of PP-, PA- and AA-butterflies with Z a = 75, Z b = 15, left panels and Z a = 15, Z b = 75, 
right panels. Initial conditions are given by Eqs. (34a), (35) and (36). In Figs. 11 - 16 time is measured in arbitrary units 



c. PA- and AF '-butterflies. Considering free evolu- 
tion from asymmetrical initial conditions we will distin- 
guish PA-butterfly, in which the leading triad has P- 
conncction and the driven triad has A-connection, from 



AP-butterfly, in which the leading triad has A-connection 
and the driven triad has P-connection. 

We found that in AP-butterfiies the energy transfer to 
the driven triad is blocked by the second of the conser- 
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vation laws (27b) and its evolution is similar to that of 
PP-butterfly. 

On the contrary, PA-butterfhes demonstrate time evo- 
lution similar to that of AA-butterfiies: compare in 
Fig. 11 panels C and V for PA-butterfiy with panels 8 
and T for AA-buttcrflies. One sees that corresponding 
plots in both panels are quantitatively the same: in the 
left panels (with Z a — 75 , Z b = 15) the energy trans- 
fer is strongly suppressed (with J ~ 0.01 in the driven 
triad), while in the right panels (with Z a = 15 , Z^ = 75) 
the energy transfer is efficient (with J ~ 10 -f- 30 in the 
driven triad). 

To complete discussion on how the energy transfer de- 
pend on the ratio Z a / Z\, we present in Fig. 13 time evo- 
lution from the same as in Fig. 11 initial conditions, but 
with equal values of the interaction coefficients, taking 
for concreteness Z a = Zf, = 15. The excitation level of 
the driven triad «7i 216 — 0.2. As expected, this level is 
larger than ~ 0.01 for Z a J Z\, = 5 and smaller than ~ 20 
for Z a /Z b = 1/5. 
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FIG. 13: Color online. Suppression of the energy transfer in 
PA-butterfly with equal values of the interaction coefficients: 
Z a — Zb = 15. Initial conditions are the same as in Fig. 11 
for PA-butterfly: Eqs. (34a), (35c) and (36). 



3. Effect of initial conditions 

As we mentioned, the time evolution and efficiency of 
the energy transfer between triads crucially depend on 
the initial conditions. To demonstrate this we present 
in Fig. 14 the evolution of PA-butterflies ( panel A) and 
AA-butterfiies (panel £>) for the same "efficient" values 
of the interaction coefficients (Z a = 15, Zf, = 75), as 
in Fig. IIP, £ and T. The only difference is in initial 
conditions, that on the first glance are very similar to 
that used above in Fig. 11. Namely, for PA-butterflies we 



simply interchange initial amplitudes of individual modes 
in the leading a-triad, replacing Bo <^4> Bo in Eq. (35c). 
For AA-butterflies we replace almost equal values of Bo 
and Bo in Eq. (34a) by different in an order of magnitude 
values (34b) with the same sum l^ol 2 + |i?o| 2 ~ 30. The 
result of this "minor" change is crucial: in PA-butterfly 
the excitation level Ji^\b decreases from ~ 20 to ~ 0.4 
and in AA-butterfly from from ~ 20 to ~ 0.1. 

In order to rationalize this effect notice that the energy 
flux from the a- to 6-triad is proportional to the level of 
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conditions (34a), that leads to the best energy transfer, 
in which, according to Fig. 11 T max J\ <m, ~ 30. 

Let us show, that this is indeed the maximal possible 
value of the level of excitation of the driven triad J\ : 2\b- 
To this goal consider the invariants (27c) of AA-butterfly, 
that can be combined as follows: 



/aa = |i?i| a | 2 + |i? 2 | a | 2 +2|S3|a| 2 +|5 1 | fe | 2 + |i?: 



2\b\ 



(41a) 



At t = Bu a = Bq, B 2 \ a = Bo (or vise versa) and 
much larger than the rest of the amplitudes. Therefore 
Iaa = Bi|a| 2 + |5 2 |a| 2 - Ji,2\a = \B \ 2 + \B \ 2 and thus 

J li2 | 6 ~ |5i| 6 | 2 + |B 2 | 6 | 2 </aa«30, AA. (41b) 

Finally notice, that in view of the restriction (40) it 
is clear, that with the second choice of the initial con- 
ditions (34b) the energy transfer should be much less 
efficient. This is exactly what one sees in Fig. 14 B. 

b. PA-butterfly. The situation is a bit different for 
PA-butterfly with P- and A-modes being individual in 
the leading a-triad. In this case combining invari- 
ants (27b) for PA-butterfly one finds new (dependent) 
invariants: 



Ipa — \B 3 l a \ 2 - \B 2 \ a 



-\B 



2\b\ 



J; 



3,2|a 



2|B 3 |6l 

" ^1,216 " 



|Bi| b | 2 (42a) 



2\B. 



3\b\ 



At t = 0, B 2 \ a = B , B 3 \ a = £> (or vise versa) and 
much larger than the rest of the amplitudes. Therefore 
Ipa = ^2,3|a(0) which is equal to -(|£? | 2 - |-So| 2 ) for the 
old initial conditions and to I PA = \Bq\ 2 — \Bq\ 2 for the 
new ones. This allows one to write: 



FIG. 14: Color online. Suppression of the energy transfer by 
initial conditions in PA- and AA butterflies. 



excitation of the common mode, (|i?3|b| 2 ) for both butter- 
flies under consideration. To estimate the upper bound 
for (|£> 3 |fr| 2 ) wc can approximate a-triad as an isolated 
one, neglecting the feedback effect of the driven triad on 
the leading one. In this approximation we can use the 
invariants (19) for isolated triad: 



il,3|a= |£l|a| 2 + |5 3 |a| 2 l. 



2.3|(i 



|S 2 |a| 2 +|S 3 |a| 2 . (39) 



a. AA-butterfly. In AA-butterfly initially the ampli- 
tude |fl 3 | | < \B 1{a \, \B 2 \ a \. Therefore 7 lj3 | a = \B 1]a (t = 

0)| 2 ee \B^\ 2 and / 2 , 3 | = \B™\ 2 and thus 



\B 3{a \ 2 <\B^\ 2 
\Bs\ a \ 2 <\<! ]2 



|B 3 | a | 2 <min{|^| 2 ,|i?W| 2 }.(40) 



This means that the efficiency of the energy transfer is 
determined by the smallest initial magnitudes of the in- 
dividual modes; fixing sum of their squares one has the 
most efficient transfer at (almost) equal initial magni- 
tudes. This was realized in the first version of the initial 



/pa = J 2 , 3 |a(0)«±26.7, PA. 



(42b) 



Unfortunately, in this case invariant I rA is not positive 
definite and therefore one cannot get rigorous restriction 
on Ji,2\b> similar to (41b). Nevertheless, more detailed 
analysis allows us to think that the upper bound for J\^\b 
for PA-butterfly should be the same as for AA-butterfly, 
i.e. about 30 for our initial conditions. 

In order to clarify the dependence of the energy trans- 
fer on initial conditions we should estimate the amplitude 
of common mode similarly to the case of AA-butterfly, 
neglecting the excitation of the driven triad. In this case 
invariants (27b) can be written as: 



3\a 



l|ol 



I^Slal 2 h 



2\a 



\B 



l\a\ 



\B 2 



(43) 



At initial moment of time the amplitude |-Bi| a 

\B 2 \ a \, \B 3la \. Therefore J 1)3 | a ~ |B 3 |„(i - 0)| 2 

and I 3i2 \ a = \B^\ 2 - \B^\ 2 . Thus the first of Eqs. (43) 
gives 



|d(0)|2 
\B 3 \ a I 



\B lla \ 2 = \B 3lb \ 2 <\B^\ 



(44) 



while the second one leads to the trivial restriction 
|-Bi|a| 2 > — l-^3|o l 2 ' tnat ^ s satisfied automatically. The 
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FIG. 15: Color online. Time evolution of AA-butterfly with different level of excitation in the driven triad, governed by 
parameter C in Eq. (36). 



conclusion is that the efficiency of the energy transfer is 
determined by the initial magnitudes of the individual 
A-mode. This was realized in the first version of the ini- 
tial conditions in which we assigned more energy to B^\ a 
mode; that leads to the best energy transfer. 

The next question is how to rationalize why the energy 
transfer into the driven triad can exceed 50%. Intuitively 
the answer is rather obvious: we understood already that 
the energy transfer is much more effective, if the accept- 
ing triad has larger interaction coefficient Z. Therefore 
in the considered case, when Z\, Z a , during long evo- 
lution with various values of the triad phase (that deter- 
mines the direction of the energy flux) and with similar 
value of the triad excitations, the energy flux from a- to 
6-triad is more favorable than the flux in the opposite di- 
rection. This leads to the asymmetry of the mean energy 
content between triads in favor of the 6-triad with larger 
value of the interaction coefficient: Zf, ^> Z a . 

c. What depends on the excitation level of the driven 
triad? Up to now we have considered initial conditions 



in the leading triad, just by mentioning that the driven 
triads have much smaller values of the initial excitations. 
The reason for this neglect is that the time evolution 
is insensitive to the initial excitations in the driven tri- 
ads, their level just have to be very small. To demon- 
strate this we compare the time evolution of AA-triad 
with Z a = 15, Zi, = 75 from the initial conditions (34a), 
(35b) and with "standard" level of initial excitation in 
6-triad [C = 1 in Eqs. (36)] with that for the much 
smaller initial energy contents in the driven triad, see 
Fig. 15. The long-time evolution is practically indepen- 
dent of C as shown in Fig. 15 A The C-dependence is 
visible only on the initial stages of the evolution: panel-S 
with C = 1, -C with C = 0.1 and -V with C = 0.01. One 
sees that the position of the first maximum for C = 1 
is at t max « 0.03, for C = 0.1 at i max w 0.04, while 
for C = 0.01 at £ max f=a 0.05. This dependence is quite 
understandable: according to Eqs. (25b) initial small per- 
turbations in the driven triad grow exponentially in time 
due to the parametric instability of the common 
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mode with respect of decay into small individual modes: 
Bi\b , B 2 \b oc exp(Zh|S 3 |f ) |t) modes. Therefore 

imax OC -InZftC, (45) 

as observed. 



C. Energy-junctions in the triple-triad clusters 

In this Subsection we study free evolution of the triple- 
triad clusters, that allow to shed light on the energy 
transfer from one triad to two other triads. To this aim 
we consider here free evolution of various triple-triad clus- 
ters, in which the leading a-triad will be initially highly 
exited, while the level of excitation of two other b- and 
c-triads will be relatively small. One topological option 
for effective energy flux from a- to b- and c-triad, consid- 
ered in the next Sec. V C 1 is PAA-star. One more option 
is the triple-triangle configuration, in which energy goes 
from one triad to two other, PP-connected triads. This 
structure is very rare and will not be discussed here. Next 
option, considered in Sec. V C 2, is AP-PA chain with the 
leading a-triad in the middle. 



1. Triple-star junction 

An effective component in the energy transfer through 
big clusters is the PAA-star. PAA-star can accept en- 
ergy via A-mode of a-triad and transfer it to two other 
b- and c-triads, via their A-modes, B 3 \ b and _B 3 | C , con- 
nected to the (same) passive B^ a mode of the a-triad 
(see Fig. 18, A). In order to clarify how an additional 
c-triad affects the energy in PA-butterfly studied above, 
we present in Fig. 16 an evolution of PAA-star with the 
choice of the interaction coefficients similar to that of 
above and with "energy-transfer effective" initial condi- 
tions, given by Eq. (35c) in a-triad, Eq. (36) in o-triad. 
For c-triad we took initial conditions similar to that for 
o-triad Eq. (36), but with the complex conjugated Bi o 
and i?2,o- Hence, the initial conditions are as follows: 

B Ma (0) = B 3]b (0) = B 3]c (0) = £? 3 ,o , f46a) 
B 2 \ a (0)=B Q , B 3]a (0) = B ; 

B 1]b (0) = B h0 , B 2 | 6 (0) = B 2fl ; (46b) 

5i|c(0) = Bio , B 2 , c (0) = Bio . (46c) 

Quite expectedly, when the leading triad has much 
smaller interaction coefficient than the driven ones (see 
Fig. 16 A with Z a <C Z b = Z c = 75), the energy transfer 
to both driven triads is fully effective. The new clement 
with respect to the butterfly case is the energy oscillation 
between driven triads. 

Panel B in Fig. 16 demonstrates an evolution of the 
PAA-star, in which driven c-triad has the interaction co- 
efficient Z c = 3, i.e. five times smaller than in the leading 



triad. As expected, the energy transfer to very slow c- 
triad is strongly suppressed, below the level J\^\c ^ 1- 

The same level of excitation, J\^\b ~ J\.2\c ^ 1 demon- 
strates the PAA-star with all equal interaction coeffi- 
cients, see Fig. 16 C. The energy alternations between 
equivalent triads is even more pronounced, than in case 
Z a <C Z b = Z c = 15, that demonstrates a more stochas- 
tic behavior. 

The comparison of Fig. 16 B and T> shows that the 
decrease in Zb (from 75 to 15) suppresses energy transfer 
to the c-triad even further up to level 0.015 <C 1. 

The overall conclusion is that qualitatively the dif- 
ference in the interaction coefficients affects the energy 
transfer in triple stars similarly to that in butterflies, but 
on the quantitative level the interplay of two driven triads 
cannot be neglected. 

Another expected conclusion is that the triple clus- 
ters demonstrate "more random" evolution than but- 
terflies. However, the level of randomness strongly de- 
pends on the initial conditions. One sees this from three- 
dimensional parametric representation of PAA-star tra- 
jectories in Fig. 17 with different choice of the interaction 
coefficients in the upper and lower panels. Left panels 
shows randomization of trajectories, starting from com- 
plex initial conditions (46), right panels - much more 
regular behavior of the trajectories, starting from simi- 
lar, but real initial conditions (in which complex number 
in Eq. (46) are replaced by their absolute values). The 
qualitative explanations of this difference is very simple. 
Equations of motion for triple stars (29) are real. There- 
fore the mode amplitudes B , that start from the real ini- 
tial conditions, remain real during the evolution and the 
dimensionality of the phase space in that case is smaller 
than that during evolution from general (complex) ini- 
tial conditions. Note also that the real trajectories, al- 
though close to periodical ones, have small random com- 
ponents that lead to the finite "width" of the attractors 
in Eq. (46), right panels. 



2. AP-PA-chain energy junction 

Another example of an effective energy junction is the 
triple AP-PA-chains, in which energy can be accepted 
via A-mode of the middle a-triad and transferred to other 
two b- and c-triads via their A-modes, B 3 \ b and i?3| c , con- 
nected to different P-modes of the a-triad from the left 
and the right: B 3 \ b = B x \ a and B 3 \ c = B 2 \ a (see Fig. 18, 
B) . A choice of the initial conditions that guarantees ef- 
ficiency of the energy transfer is similar to Eqs. (46) with 
the differences that are dictated by different types of the 
connections: 

Bi| (0)=B 3 | 6 (0)=i?3,o, (4?a) 
B 2 \a(0) = B 3 |c(0) = Bio, B 3la (0) = Bo; 

5i|b(0) - Bi l0 , B 2lb (0) = S 2j0 ; (47b) 
B 1]c (0) = B* h0 , fl 2]c (0) = 5 2 *, . (47c) 
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FIG. 16: Color online. Efficiency of PAA-star energy junction with different Zb/Z a and Z c /Z a [initial conditions (46)]. 



In Fig. 19 we show the time evolution of the triple 
AP-PA-chain with the same four sets of the interaction 
coefficients, as in Fig. 16 for APP-stars. 

Comparing Fig. 19 with Fig. 16 panel by panel one con- 
cludes, that the triple AP-PA-chain, in which the driven 
triads are connected to different modes of the leading 
triad, demonstrates much more efficient energy transfer, 
than the triple PAA-star, in which the driven triads are 
connected to the same mode of the leading triad. An- 
other difference, is that the time evolution in the triple 
AP-PA-chain is much more intermittent than in the triple 
PAA-star. 



VI. FINITE-DIMENSIONAL WAVE 
TURBULENCE IN THE LONG-CHAIN 
CLUSTERS 



In this Section we study stationary energy distribu- 
tions of resonant waves and energy exchange between 
triads in the long-chain clusters in the regime of finite- 
dimensional wave turbulence with constant energy flux 
and during free evolution. 
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FIG. 17: Color online. Three dimensional parametric representation of trajectories of PAA-star with Z a = 3, Zb = 15, Z c — 75, 
panels A and and with Z a — 3, Zt = 15, Z c = 75, panels C and X>. Left panels: Complex initial conditions (46); Right 
panels: real initial conditions in which complex number in Eq. (46) are replaced by their absolute values. 
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FIG. 18: Triple-star (.4) and triple-chain (B) energy junc- 
tions. Initial energy is localized in the leading a-triad and 
goes during free evolution in two driven b- and c-triads. 



resonant triads (with N = 12, 16, 20 and 24). To do this 
numerically one should model the energy pumping in the 
leading a-triad and energy dissipation in the driven N th - 
triad. 

The simplest and reasonable way to model the energy 
dissipation is quite obvious: one should introduce linear 
damping terms into the equations of motion for individ- 
ual modes of the last driven triad. In our case these are 
P-modes and instead of Eqs. (25c) we suggest the follow- 
ing two equations for them: 



B 



1\N 



B 2 \n — ZnB*, n B 3 \ 



N 



1i\nBi\n i 

72\nB2\N , 



(48) 



A. Pumping and damping in chain clusters 

With this subsection we begin to consider the station- 
ary dynamics of the A^-chain clusters, consisting of many 



with some phenomenological damping frequencies 7i|at 
and 72 1 at that should be compared with the characteristic 
interaction frequency in these equations 



^intlN = Z 



N 



IB 



31AM 



(49) 
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FIG. 19: Color online. Efficiency of AP-PA-chain energy junction with different Zb/Z a and Z c /Z a [initial conditions (47)]. 



Then one can distinguish the cases of symmetrical damp- 
ing: 7i|iv = 72|iVj asymmetrical damping: ^ 7 2 |jv, 
but 7i|jv — 72 1 at and strong damping asymmetry: 7i|jv 3> 
7 2 1 TV- I n addition there are cases of weak damping, when 
\Jl1\Nl2\N ^ ^int|N 5 intermediate and strong damping, 
when ^7i|jv72|at ~ ^int|N and ^n\Nl2\N > ^int|N- To 
start with in this paper we will restrict ourselves by the 
simplest case of weak symmetrical damping. 

In all cases we will include energy pumping term 
in the equation of motion for the A-mode of the first 
(leading) triad, i.e. in the equation for B^\ a . Now it 



reads: 

B m =-Z 1 B 1 \ 1 B 2 \ a +P 3 \ 1 . (50) 

How to mimic energy pumping is a much more delicate 
question. Notice, that in the problems of developed wave 
or hydrodynamic turbulence researches usually fill free to 
model the energy source and sink in the simplest possible 
manner [27]. For example, one introduces a random force 
with Gaussian statistics (that has some justification) and 
(5-corrclatcd in time, which is usually far from reality, 
see e.g. [28]. The rationale is that for very high levels 
of turbulence excitation the incrtial interval of scales is 
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large enough such that one expects universal statistics of 
turbulence, independent of characteristics of the energy 
pumping [27]. In our case of relatively low level of exci- 
tation, when the number of excited degrees of freedom is 
not so huge, the universality is, generally speaking, ques- 
tionable: statistics of the system dynamics can depend 
on the way the system is excited. One can imagine few 
very different versions of the pumping term. 

The first version mimics an instability of the B 3 ^ mode 
with an inverse growth time frail) : 

^3|i = ^3|i^3|i i (51a) 

Next one can consider a periodic driving force: 

P 3 |i = F u cos(wt) , (51b) 

that can mimic the effect of an "external" triad, con- 
nected to i? 3 |! mode. In this case w is the leading fre- 
quency of energy exchange between modes of the external 
triad. 

Also one can take a random white Gaussian noise 
Pan = /(*), (f(t)f(t'))=f 2 S(t-t'), (51c) 

that mimics the influence of numerous non-resonant tri- 
ads. 

In this paper we will utilize a bit more realistic driving 
force P3\ a (t) [periodic, but not oc cos(wf), as in Eq. (51b)] 
that is produced by the equations of motion for an iso- 
lated triad. This mimics the effect of an "external" triad, 
connected to B 3 \ a mode in the approximation, when the 
feedback effect is neglected. 

B. Condition of stationarity 

Consider first the dynamical invariants for isolated N- 
chain clusters (without pumping and damping terms). 
Combining invariants (19), one has in this case a (depen- 
dent) invariant for isolated triad in the form 

'|i = |5 3 |i| 2 + i(|Bi|i| 2 + |5 2 |i| 2 ) ■ ( 52a ) 

Similarly, a combination of invariants (27b) yields the 
invariant for PA-butterfly: 

I\i,2 = |£ 3 |i| 2 + \B 3 \ 2 \ 2 + l{\B ll2 \ 2 + \B 2{2 \ 2 ) . (52b) 

And finally, from invariants (31c) one derives new (de- 
pendent) invariant for PA-PA-chain: 

I\i,2,3 = \B3\1\ 2 + \B 3 \ 2 \ 2 + \B 3 \ 3 \ 2 + -(|Bi| 3 | 2 + |-B 2 | 3 | 2 ) ■ 

(52c) 

Similarly, for isolated iV-chain cluster one gets: 
N 1 

/|i,2,..jv = E l B % ! 2 + 2 ( |Sl l w|2 + |jB2 I w|2 ) • (52d) 
J'=l 



All these invariants include B^-modc, affected by the 
pumping, and ^ ( | jv 1 2 + 1 1 jv 1 2 ) terms, affected by the 
damping. Therefore under the stationary conditions the 
following equality must hold: 

(P3\iB* 3 \i) = \(li\N(\B llN \ 2 )+ l2lN (\B 2lN \ 2 )) . (53a) 

In particular for the pumping (51a) one has a simple 
relationship: 

73|i<|5 3 |i| 2 ) = 1{ii\n(\B mn \ 2 )+ 121n (\B 21n \ 2 )) • (53b) 

Equations (53) can be considered as the (nessesary) con- 
dition of the energy balance: the left hand side describes 
the energy pumping, while the right had side - the energy 
damping. Each of these quantities can be equated (under 
the stationary conditions) to the energy flux through the 
clusters. 



C. Universality of statistics of finite-dimensional 
wave turbulence in the long-chain clusters 

1. Modeling the parameters of the forced chain clusters 

We present here our choice of the modeling parameters 
for the forced PA-PA-. . .-PA chain clusters. The cluster 
consists of N triads and is constructed by connecting the 
first P-mode of (n — l)'s triad, £i| n _i, with the A-mode 
of the next, n's triad, B 3 \ n . The interaction coefficients 
Z n were chosen in geometric progression, 

Z n = AZ„_! = Z x \ n - X . (54) 

The chain is forced by an independent freely evolving 
triad with Zq = 0.125 related to Z\ as Z\ = Zq\. The 
amplitude B\\n of the forcing triad is added to the equa- 
tion of i?3 1 1 of the first triad in the chain 200 times during 
the period 2/ of the forcing triad. We made special ef- 
forts to ensure the exact periodicity of the forcing. 

The linear dissipation with the damping frequency 7 
was added to the equations of two P-modes of the last 
triad, B^ N and B 2 \ N . In most calculations 7 = 10~ 4 was 
used. 

We have studied the mean square amplitudes of "free" 
modes, (|i?2|n| 2 ) and "connected" modes (|i? 3 |„| 2 ) = 
(|-Bi|„_i| 2 ) in the chains with different number of tri- 
ads, N = 16, 20, 24, for the set of the frequency scaling 
parameter A = 1, y/2, 2, 2v / 2 and using different initial 
conditions in the pumping triad. 

We found that initial conditions do not effect the sta- 
tionary statistics in the cluster. For concreteness, for all 
presented simulations we chose for the pumping triad 

B 1{0 = 0.05 + 0.02 i, (55) 
B 2l0 = 0.02 + 0.05 i, (56) 
B 3 \ = 0.9- 0.93 i; 
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and for all TV triads of the chain 

B 1]n = 0.05 + 0.02 i, (57) 
B 2]n = 0.02 + 0.05 i, (58) 

£>3|ra — -Bl|n-1 • 

The time of the chain evolution required to reach sta- 
tionarity varied from about 1000 Tf of the pumping triad 
for A = 1 to about 50 Tf for A = 2^/2. For times later 
than 1000 Tf we calculated mean values in time, aver- 
aging values of (\Bj\ n \ 2 ) over 3001/ . In the most fast 
converging case of A = 2V2 we averaged over 50T/ . In 
all cases the evolution was followed for much longer times 
and we verified that a particular choice of the averaging 
interval does not affect results. 

Our results are presented in Figs. 20 as logarithm of 
(\B 2 \ n \ 2 ) and (|B 3 |„| 2 ), normalized by its value in the first 
triad vs triad number n normalized by TV. To force all 
lines to start from the same point and not from 1/TV, we 
shifted triad numbers by 1: n — > n — 1 for presentation. 
Such a shift leads to better collapse of the data. 

2. N-, X- and 7 -independence of the forced-chain statistics 

Figure 20 ^4 presents normalized mean-square ampli- 
tudes (\B 2 \ n \ 2 ) an d (|-B 3 | n | 2 ) for the chain with TV = 16 
and a set of the scaling parameters A = 1, y/2, 2 and 
A = 2y/2. One sees that the lines for different A arc quite 
close, quickly converging for larger A. We conclude that 
the distributions of (|i?2|n| 2 ) and (|i? 3 |„| 2 ) are almost A- 
independent at least for A > 2. 

Figure 20 B represents distributions of (|-E>2|n| 2 ) and 
(|i?3|„| 2 ) for the chains with A = \/2 and different number 
of triads: TV = 12, 16, 20 and TV = 24. One sees that the 
lines for different TV practically coincide. The conclusion 
is that the distributions of {\B 2 \ n \ 2 ) and (|B 3 |„| 2 ) are TV- 
independent in the limit of large TV. 

Panel C in Fig. 20 shows (|S 2 |„| 2 ) and (|-B 3 |„| 2 ) for the 

chains with A = s/2, TV = 16, and different, but small 
values of 7 = 10~ 4 , 5 x 10~ 4 , 10~ 3 , 5 x 10~ 3 . Collapse 
of these plots is obvious. One can ask, what happens if 
one increases further the damping parameter? Increasing 
7, we found that in this case the individual amplitudes 
in the last triad become very small. Its ability to absorb 
the energy flux from the previous triad diminishes to the 
stage at which it can not anymore serve as an energy 
sink. This leads to accumulation of the energy in the 
N-l'th triad and later in N-2nd triad etc, thus creating a 
"bottleneck effect". 

The reason for this effect is that the energy pump- 
ing of Si|jv and B 2 \n modes is oc Ke[B^ N Bi^B 2 ^~\, 
while their energy damping (in the simple case 7x ijv = 
72, jv = 7) is proportional to 7(|Bi|at| 2 + |-B 2 |aH 2 )- In 
other words, both fluxes are quadratic in the amplitudes 
Siijv and i? 2 |7V ■ It means that for large enough 7 the 



energy damping will be larger than the energy pumping, 
whatever the values of B^ N and B 2 \ N are. In this case 
these amplitudes will decrease in time exponentially. In 
the other case, when 7 is small enough, the energy damp- 
ing will be smaller than the energy pumping and ampli- 
tudes B^ N and B 2 \ N will exponentially grow. In other 
words, there is a critical value of 7, above which the last 
triad is not exited at all. A more detailed analysis shows, 
that for 7i|jv 7^ 72. n the critical parameter is ^/7i|at 72, n- 



3. Forced- vs free- evolution statistics in the chain clusters 

One more surprise is that in the conditions of the 
energy-flux equilibrium in the scale-invariant situations 
(in our case, Z n oc A") we did not observe the expected 
power-like behavior, when (|_Bj|„| 2 ) oc A"^. Moreover, 
the ratio (\B 2 \ n \ 2 ) f (\B^ n \ 2 ) depends on n and the con- 
nected modes { | -B3 1 n | 2 ) weakly depend on n, as one ex- 
pects in the thermodynamic equilibrium, after a long free 
evolution. 

That is why we decided to compare forced distribu- 
tions of (|i?2|n| 2 ) and (l-B^I 2 ) with the distributions of 
(l-^2|n| 2 ) and (|i? 3 |„| 2 ) after long free evolutions (without 
forcing and damping) in the same cluster. For concrete- 
ness we took cluster with A = \/2 and TV = 16 after forced 
evolution during time t = 1000T/. To get the values of 
the forced amplitudes we averaged over last 500T/. Then 
we took final complex amplitudes B 2 \ n and _B 3 | T1 as initial 
conditions for the free evolution during another period 
t = 10007/ , and again averaged over last 500T/. Both 
distributions of (|7?2|n| 2 ) and (|7? 3 | n | 2 ) are presented in 
the Fig. 20 T>. As one sees these distributions practicably 
coincide. 

In other words, the distributions with constant energy 
flux and the distributions with zero energy flux are al- 
most identical. To show, why this happens we presented 
in Fig. 21 the time evolution of the instant energy flux 
via n-th triad 

e n (t) = 2Z n Re[B lln (t)B 2ln (t)B* 3 \ n {t)} , (59) 

normalized by its mean value e n = (e n (t)}. One sees 
cnourmous fluctuations of £„ (t) with peaks, reaching as 
6000 e rl . The computed mean value of the flux fluctuation 
6e„ = \J (e 2 (tf) is about 100 times larger than e n . This 
means that there is a strong energy exchange between 
triads in the cluster, that essentially exceeds minor mean 
flux. Therefore one can switch off the mean energy flux 
without almost any effect on the cluster statistics. Wc 
think, that this general statement is valid for all kind on 
clusters, consisting of connected triads. 

The only difference between forced and free long-time 
evolution of the clusters is the restriction, that follows 
from the conservations laws, that are satisfied exactly in 
free evolutions and approximately in the forced case. 
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FIG. 20: Comparison of the plots of In (|_Bj| n | 2 ) vs nj [N — 1) for different A, panel A, for different N, panel B and for different 
7, panel C. Panel V compares distributions for the free and forced evolutions. 



VII. CONCLUSIONS 

A. Main points of understanding in 
finite-dimensional wave turbulence 

• In the first part of the paper we studied the structure 
of finite clusters of resonant triads using the example of 
atmospheric planetary waves and showed that: 

- In the physically relevant domain of atmospheric 
planetary waves (m , £ < 1000, when the mode-scale 



is larger than the height of the Earth atmosphere) we 
have found and described the topology of all the clusters 
formed by resonantly interacting planetary modes. The 
set od clusters contains isolated triads and subsets of 2-, 
3-, . . . , 16 and 3691 connected triads, with 2- 3-, . . . , 
9-mode and (maximum) 10-mode connections; 

- Analyzing the integrals of motion we suggested a 
classification i) of triad modes into two types - active (A) 
and passive (P), and ii) of connection types between tri- 
ads - AA, AP and PP. We showed explicitly that through 
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FIG. 21: The time evolution of the normalized flux over triad 12 (A) and over triad 16 (£>) during one period of the forcing 
triad for X = 2 for N = 20. 



the AA-connection the energy can flow in both directions, 
through the AP-connection only from P-connected triad 
to the A-connected one, but not vice-versa. The PP- 
conncctions are almost impenetrable for energy in both 
directions. Therefore from the viewpoint of the energy 
transfer, large clusters can be subdivided into smaller 
ones (by cutting PP-connections); 

- We introduced a notion of PP-irreducible sub- 
clusters that cannot be further divided by cutting PP- 
connections and studied their statistics for planetary 
waves in the spectral domain m , t < 1000. There are 
3005 triads, 400 PA- and AA-butterflies, 143 triple-triad 
PP-irrcducible sub-clusters, etc., and only two 130-triad 
ones in this domain; 

- To first approximation, almost all triads in the 
methodologically significant domain, Table. 1, can be 
considered as completely or almost separated from the 
rest of the atmospheric planetary waves and therefore 
the energy oscillations between them can really lead to 
intra-seasonal oscillations in Earth's atmosphere as sug- 
gested in [6]. 

• In the second part of the paper we presented a gen- 
eral analysis of the energy transfer in the PP-irreducible 
sub-clusters. Studying free evolution from asymmetrical 
initial conditions, when almost all the initial energy is 
localized in one (leading) triad we showed: 

- the energy flux through the AA- and PA-connections 
are qualitatively the same and crucially depend on the 
initial conditions in the leading triad and on the ratio 
of the interaction coefficients in the leading and driven 
triads; 

- the energy flux from the leading to two driven triads 
in triple-triad clusters depends on the type of energy- 
junction (triple AAA-, PAA-stars or PA-AP chain) and 
on the relations between three interaction coefficients. 

• We also studied forced stationary energy transfer in 
a long chain consisting of N PA-connected triads with 



the interaction coefficients Z n cx A™, 1 < n < N, with 
forcing in the first triad and damping in the last one. 
We showed that stationary energy distributions between 
triads are universal in the following sense: 

For A > 1 the distributions arc almost A- 
independent; For A > 2 the distributions practically col- 
lapse on each other. 

the distributions are practically independent of the 
damping parameter 7 for 7 < 7th- For 7 > 7th the 
stationary state does not exist. 

Distributions for the forced case and for the free 
evolution (with the forced state as initial condition) co- 
incide due to extremely large fluctuations of the energy 
flux, exceeding its mean value in orders of magnitude. 

• Our analysis of finite-dimensional wave turbulence 
is based on the simple form (18) of the basic equations 
and does not exploit the explicit form of the interaction 
coefficients Z which vary substantially for different wave 
systems. This means that our results can be used di- 
rectly for arbitrary 3-wave resonant systems governed by 
Eqs. (18), e.g. drift waves, gravity-capillary waves, etc. 



B. Remaining unclear issues in finite-dimensional 
and mesoscopic regimes of weak wave turbulence 

• Having in mind the extremely rich variability of the 
structure of concrete finite clusters of resonance triads 
and the various important aspects of the cluster dynamics 
and statistics (some of them were not mentioned at all) , 
this paper opens more questions than answers. Among 
them: 

- What might be the peculiarities of free evolutions of 
finite clusters from other initial conditions, more general 
than those studied in this paper, especially for particular 
clusters which may be important in some applications? 
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How energy goes through more complicated junc- 
tions, than triple stars, studied in this paper, such as 
stars with four and five and more triads, etc.? 

- How integrability of clusters (with special choice of 
the ratios of interaction coefficients and/or initial condi- 
tions [25, 26]) or closeness to the integrable cases affect 
the energy transfer and the statistics of the mode ampli- 
tudes? 

What is the dimensionality of the cluster trajec- 
tories, how does it depends on the details of the initial 
conditions and/or the ratios of interaction coefficients for 
various cluster structure? 

- How different-mode, different-time correlation func- 
tions depend on the mode position in a cluster and on 
the time difference between them? 

- What is the values of the flatness (ratio of the forth- 
order correlation functions to the square of the second- 
order ones) and how does it depend on the cluster struc- 
ture, initial conditions, etc.? 

- To what extent can the statistics of the amplitudes 
of individual modes be considered close to Gaussian at 
least for very large (but finite) clusters? 

- How does this closeness (if it exists) depend on the 
mode position in a cluster and on initial conditions? 

- If this closeness exists (we believe that it does) , how 
can one formulate appropriate closures that will lead to 
an analytical statistical description of the finite cluster 
behavior? 

• A few of the most important questions, at least from 
the theoretical viewpoint, are: 

What is the principal difference between finite 
size cluster behavior for three- wave resonances, discussed 
here, and that for four-wave resonance systems? 

How mode statistics and possible applicability of 
the closure procedures changes, if one accounts for small 
damping of the mode energy and external random noise, 
that can mimic interaction of a cluster with the "rest of 
the world"? 

- How all the features of a cluster behavior (both for 



those studied and those left open) get modified if one 
accounts for quasi-resonances that can become crucially 
important with increasing the level of the modes excita- 
tion? 

How to describe the statistical behavior of finite- 
size dynamical systems (in the case of three- and four- 
wave interactions) with further increases in the system 
size or in the level of system excitation? How this behav- 
ior approachs (step by step, probably via an intermediate 
kind of behavior) the limit of infinite system with quasi- 
Gaussian statistics, that can be successfully described 
with the help of wave-kinetic equations? 



C. The road ahead 

Our feeling is that all these (and many similar and 
different, but related) questions are the subject of a new 
fields in nonlinear wave physics, finite- dimensional wave 
turbulence and mesoscopic wave turbulence. This subject 
offers very interesting issues both from the physical and 
the methodological viewpoints, with possible important 
applications in numerous physical examples. 
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